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ABSTRACT 

The ever increasing level of precision achieved by present and future radial- velocity instruments is opening the way to discovering 
very low-mass, long-period planets (e.g. solar-system analogs). These systems will be detectable as low-amplitude signals in 
radial- velocity (RV). However, an important obstacle to their detection may be the existence of stellar magnetic cycles on 
similar timescales. Here we present the results of a long-term program to simultaneously measure radial-velocities and stellar- 
activity indicators (Call, H a , Hel) for a sample of stars with known activity cycles. Our results suggest that all these stellar 
activity indexes can be used to trace the stellar magnetic cycle in solar-type stars. Likewise, we find clear indications that 
different parameters of the HARPS cross-correlation function (BIS, FWHM, and contrast) are also sensitive to activity level 
variations. Finally, we show that, although in a few cases slight correlations or anti-correlations between radial-velocity and the 
activity level of the star exist, their origin is still not clear. We can, however, conclude that for our targets (early-K dwarfs) 
we do not find evidence of any radial-velocity variations induced by variations of the stellar magnetic cycle with amplitudes 
significantly above ~1 m/s. 

Key words, planetary systems - Stars: activity - Stars: fundamental parameters - Techniques: spectroscopic - Techniques: radial 
velocities - starspots 



1. Introduction 

Following the discovery of the first extrasolar p lanet or- 
biting a solar type star ( Mavor fc Quelozl 1995), a mul- 
titude of other planetary systems ha ve been announced 
(for a review see Udrv &: Santosl 2007 ). At first, only very 
short-period companions were found, something that was 
quite unexpected by the theories of giant planet forma- 
tion. However, the impressive increase in precision of the 
current radial- velocity (RV) "machines" (the majority of 
the established exoplanets have been discovered by radial- 
velocity surveys) and the long baseline of the measure- 
ments have already brought to light both long-period plan- 
ets, much more similar to the Solar System giants (e.g. 



Wright et al. 


2008) 


Mavor et al. 


2009). 



* Based on observations collected at the La Silla Parana 
Observatory, ESO (Chile) with the HARPS spectrograph at 
the 3.6-m telescope (ESO runs ID 072.C-0096 073.D-0038 
074.D-0131 075.D-0194 076.D-0130 078.D-0071 079.D-0075 
080.D-0086 081.D-0065). Tables 5 to 8, with the data used for 
Figs. 2, 3, and 8, are available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http: / / cdsweb.u-strasbg.fr / cgi-bin / qcat? J / A+ Aj\ 



Although extremely efficient, the RV technique is an 
indirect method, and it does not allow direct detection 
of a planet, but only measures the (gravitational) stellar 
wobble induced by the planetary companion orbiting its 
sun. As a consequence, one of the problems is that one 
has to assure that the radial-velocity variations observed 
are not being caused by some other mechanism unrelated 
to the presence of a low- ma ss companion. Phen omena 
such as stellar pulsation (e.g. Bouchv et al.l 2004 ), inho- 
mogeneous convec t ion, or spots ( Saar fc Donahue! Il997 ; 
Santos et al. 120001 IPaulson et al.l l2002h can prevent us 
from finding planets (if the perturbation is stronger than 
the orbital radial- velocity variation); but perhaps more 
important, they might g ive us false candidates if they pro- 
duce a periodic signal dQueloz et al.l I2001 ; Bonfils et al 



l2007t iHuelamo et al.l I2008T I. In other words, the radial- 
velocity technique is sensitive not only to the motion of a 
star around the center of mass of a star/planet system, but 
also to possible variations in the structure of the stellar 
surface. 
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Table 1. Basic data for our stars and observation log. 



Star V B— V Sp. Type P rot [days] First observation Last observation # nights <S/N> 



HD4628 


5.74 


0.890 


K2V 


39 


2003-10-26 


2008-09-05 


36 


34 


HD16160 


5.79 


0.918 


K3V 


48 


2003-10-26 


2008-09-05 


43 


28 


HD26965A 


4.43 


0.820 


K1V 


43 


2003-10-26 


2008-09-05 


38 


43 


HD32147 


6.22 


1.049 


K3V 


47 


2003-10-26 


2008-09-05 


35 


24 


HD152391 


6.65 


0.749 


G8V 


11 


2003-06-18 


2008-09-11 


38 


36 


HD 160346 


6.53 


0.959 


K3V 


37 


2004-02-10 


2008-09-11 


31 


29 


HD216385 


5.16 


0.487 


F7IV 


7 


2003-10-26 


2008-09-11 


29 


56 


HD219834A 


5.20 


0.787 


G6/G8IV 


42 


2003-10-26 


2008-07-09 


25 


48 



This matter is getting more important every day. One 
clear output of the increase in precision of the current 
radial-velocity surveys is the ability to find low-mass 
and long-period planets, that are able to induce long- 
period but low-amplitude radial-velocity signals. These 
expected Jupiter-analogs might be easily fo und e.g. with 
instruments like HARPS (jPepe et al.ll2002h and with fu- 
ture instrumentation like th e projected ESQ inst ruments 
ESPRESSC0 and CODEX (IPasauini et al.ll2008h . But if 
in these cases stellar rotational effects are not potentially 
dangerous (rotational periods are much shorter than the 
long-period Jupiter orbit), another almost quantitatively 
unknown effect might threaten these detections: stellar ac- 
tivity cycles. 

It is now well known that many solar type stars hay e 
important magnetic activity cycles (jBaliunas et al.lll995l ). 
similar to the solar 22-year period (11+11 period spot 
cycle) . These cycles have timescales comparable to the or- 
bital periods of real Jupiter like planets. In this sense, RV 
measurements of the Sun during a 5-year period have re- 
vealed that "our" star might be s table up to a precision of 
a few m/s (jMcMillan et al.lll993l ). However, no tests have 
been done regarding other solar-type stars. Some effect of 
the activity level on the radial- velocity might indeed be ex- 
pected. It is known that magnetic fields are able to change 
the convection patterns (e.g. inhibiti ng convection) , thus 
changing line bisectors and line shifts (lDravinslll982l). But 
with the exception of a very few cases (jKurster et al.| [2003. 
for Barnard's star, an M dwarf), it is not really known or 
accurately explored, to which level the line asymmetries 
induced by changes in the convection patterns througout 
the stellar magnetic cycles can influence the measurement 
of precise radial velocities. 

In this paper we present a study of the relation be- 
tween long-term chromospheric activity variations, RV, 
and different line profile indicators. In Sect. [5] we present 
our sample and observations, and in Sect. [3] we derive pre- 
cise atmospheric parameters and masses for our stars. In 
Sect.|4] we then present the methodology used to derive 
the stellar activity level of our stars from each obtained 
spectrum. The analysis of the data and the presentation 
of the results is done in Sect.[SJ We conclude in Sect. El 



http:/ /espresso. astro. up.pt 



2. Sample and observations 

The choice of our targets was based on the sample 
of stars followed within the Mount Wilson project to 
study the magnetic cycle variations in near by FGK stars 
(jVaughan et al. 19781 : iBaliunas et al. 19951 ). In particu- 
lar, we based our c hoice on the sample presented in 
Baliunas et al. (1995). From the stars and results pre- 
sented by these authors, we first took those showing clear 
activity cycles. Then, only objects at a declination south of 
+10 degrees were considered, in order to be able to follow 
them from the La Silla observatory (the site of HARPS) 
for several months each year. One star that presented a 
particularly stable activity level was also chosen as the 
standard (HD 216385). 

In Table[T] we list the stars in our sample with their 
V magnitudes, B-V colors, and spectral typ es. The va l- 
ues were taken from the Hipparco s catalog ( ESA 1997). 
Rotational periods were taken from lBaliunas et al. r il996li . 
As we can see from the table, most stars are late-G or 
early-K dwarfs, with the exception of HD 216385 (late F) 
and HD219834A (likely a subgiant) - see also Table[2j 

To allow for a continuous follow-up of each target, 
observations were done in service mode starting from 
September 2003 (ESO period 72), just when HARPS be- 
came available. All stars were followed until September 
2008 (ESO period 81). No more follow-up was then pos- 
sible due to the end of service-mode observations at La 
Silla. For HD 152391, one extra measurement done during 
the commissioning of HARPS (June 2003) was added. 

In general, each star was observed in 5 different epochs 
during each observing season (6 months). An effort was 
made to spread the observations in time as much as possi- 
ble. This was important to allow us to average any strong 
variations due to rotational modulation effects out. We 
are interested in studying the effects of long term mag- 
netic cycle variations. 

Every measurement was done with the high-precision 
simultaneous calibration mode. Though the use of the 
Thorium- Argon simultaneous calibration may complicate 
the subtraction of scattered light over the whole image 
(specially in the blue part of the spectrum where the spec- 
tral orders are separated by only a few pixels), this mode 
was judged important since we need m/s precision in our 
measurements. 
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Table 2. Stellar parameters derived for the stars in our sample. 



Star 


T e ff 


log g spec 


Ct 


[Fe/H] 


N(Fel/Fen) 


a(Fel/Fen) 


log ghipp 


Mass 




[K] 


[cms -2 ] 


[kms -1 ] 








[cms -2 ] 


[M ] 


HD4628 


5044±52 


4.54±0.36 


0.53±0.18 


-0.32±0.11 


257/31 


0.11/0.16 


4.61 


0.72 


HD16160 


4921±103 


4.54±0.54 


0.46±0.43 


-0.16±0.19 


258/32 


0.18/0.24 


4.58 


0.69 


HD26965A 


5136±50 


4.46±0.30 


0.44±0.15 


-0.32±0.09 


258/32 


0.09/0.13 


4.42 


0.65 


HD32147 


4902±136 


4.33±0.73 


0.60±0.54 


0.21±0.21 


259/34 


0.19/0.34 


4.62 


0.79 


HD152391 


5474±25 


4.46±0.19 


1.01±0.04 


-0.04±0.06 


256/33 


0.06/0.08 


4.56 


0.92 


HD 160346 


4998±68 


4.44±0.38 


0.81±0.17 


-0.08±0.14 


258/32 


0.14/0.18 


4.63 


0.78 


HD216385 


6371±30 


4.24±0.34 


1.63±0.04 


-0.13±0.06 


230/36 


0.06/0.12 


4.06 


1.36 


HD219834A 


5530±41 


3.95±0.17 


1.06±0.04 


0.13±0.11 


257/34 


0.11/0.07 


3.96 


1.25 



Stellar oscillations can induce significant radial 
velocity signals on timescales of a few minute s. In order to 



averag e the stellar oscillation modes out (e.g. lSantos et al 
2004aj), each measurement was done with a total exposure 



time of 15 minutes. For the brightest stars in our sample, 
this implies several shorter exposures to avoid CCD satu- 
ration. 

Individual RVs were derived with the HARPS pipeline. 
The velocities were derived using the cross-correlation 
function (CCF) method. The bisector inverse slope 
(BIS) was derived using the methodology described in 
Queloz et al.l (|2000t . Other parameters of the CCF, such 
as the contrast and FWHM, were also computed. 
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Fig. 1. Comparison of two spectra of HD4628 near its 
maximum (blue continuous line) and minimum (red 
dashed line) activity level in the regions of the Ca II H line 
(top), the Hq, line (middle), and the Hei D3 line (lower 
panel). The vertical dashed lines denote the regions used 
to derive the activity indexes. 



3. Spectroscopic stellar parameters and masses 

With the combined HARPS spectra, we derived stellar 
parameters and metallicities for the sampl e stars using 
the me thodo logy and lin e - lists described in ISantos et al.l 
(2004b) and ISousa et al. ( 20081) . respectively. Equivalent 
widths of individual Fel and Fell lines were measured 
with the automatic code ARElfl A grid of H uruc 3 (|l993h 
model atmospheres was adopted together with t he 2002 
versio n of the radiative transf er code MOOG (Sneden 
19731 ). We point the readers to ISantos et ail (|2004bl ) for 
more details on the technique. In Table[2]we list the final 
derived values, together with the number of Fel and Fen 
lines used and their dispersion. For comparison, the as- 
trometric surface gr avity was also derived using Eq. 1 in 
Santos et al.l (|2004bl ). Both estimates of the surface grav- 
ity agree well within the errors. 

Stellar masses were also deri ved by int e rpolat - 



ing the theoretical is ochro nes o f Schall er et al.l (|1992T ) 



Schaerer et al. ( 1993a ). and Schaerer et al.l jl993bl). us ing 



My computed using Hipp arcos parallax es (|ESAlll997l ). a 
bolometric correction from Flower (Il996l) . and the T c ff ob- 
tained from the spectroscopy. The results are also listed 
in Tabled We estimate that the uncertainties are around 
10% because of errors in the input pa rameters and also be- 
cause of different systematic effects (jFernandes fc Santos 
20041) . 

A look at Table[2] shows that all but one of the targets 
are main sequence dwarfs. The exception is HD 219834 A, 
most likely a subgiant. 



4. Activity indexes 

From each individual spectrum, stellar activity indexes 
were derived using three different indicators: the Can H 
and K lines, the H Q line, and the He I D3 line at 5875.6 A. 
All activity indexes were derived followi ng the general pro- 
cedures described in Santos et al. I (l2000l . for the Ca 11 lines) 



and iBoisse et al.l (|2009L for H a and Hei). All the activity 
level derivations used the pipeline reduced and wavelength 
calibrated 2-dimensional spectra. Individual errors were 
derived with the number of counts in each wavelength re- 
gion. 



http:/ /www. astro. up.pt/~sousasag/ares 
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Table 3. Average ch romosphcric activity indexes derived using the Can H and K lines, with the values presented in 
Baliunas et al.l (|1995l ) presented for comparison. 



Star 


< Smw > (ours) 




<logR' HK >^ 


HD4628 


0.22 


0.230 


-4.89 


HD16160 


0.23 


0.226 


-4.88 


HD26965A 


0.19 


0.206 


-4.92 


HD32147 


0.28 


0.286 


-4.95 


HD152391 


0.39 


0.393 


-4.44 


HD 160346 


0.31 


0.300 


-4.78 


HD216385 


0.15 


0.142 


-5.00 


HD219834A 


0.16 


0.155 


-5.03 



The Can H and K "Sharps" index was computed by 
dividing the sum of the flux in two 1 A wide intervals in 
the center of each line (at 3933.66 and 3968.47 A, respec- 
tively), by the flux in two 20 A wide reference windows 
centered on 3900.0 and 4000.0 A. A weighted sum took 
the photon noise errors derived as y/~N into account (with 
N the number of counts) . Once all the spectra were ana- 
lyzed, an average Sharps value was derived for each star. 
Its value was then used to calibrate the Sharps values to 
the Mount- Wilson s cale. For this ca l ibrati on we used the 
Smw values listed in lBaliunas et al ](ll99.#l-s ee Table[3| 
Finally, using the relation in Noves et al.l ( 1984 ) , we could 
derive the values of the Ca II flux corrected for the photo- 
spheric flux, log R'hk ■ 

In Tabled we list the average values for each of our tar- 
gets (both for Smw and logi?^^), toge t her w ith the val- 
ues for the Smw listed in lBaliunas et ajj (|l995l ). As can be 



seen from the table, the values derived here agree perfectly 
with those derived by Baliunas et al. This also reflects 
that, during the 5 years of measurements, we could follow 
a significant part of the stellar magnetic cycle. Except for 
HD 152391, all the stars seem to be in the low-activity sid e 
of the Vaughan-Preston gap ( Vaughan fc Preston 1980h . 
and are thus similar to the Sun regarding their activity 
level. 

The H a index was derived by dividing the flux in the 
central 0.678 A of the H Q line (6562.808 A) by the flux 
in two reference windows between 6545.495-6556.245 and 
6575.934-6584.684 A. Finally, the He I D3 index was de- 
rived by dividing the flux in the 0.4 A central region of 
the 5875.62 A line by the flux in two 5 A wide reference 
windows in each side of the line, centered on 5869 and 
5881 A. 

While the Can H and K lines and H a line typically 
present a re-emission at the center associated with activ- 
ity phenomena, the depth of t he Hei D3 line increases 
in solar plages (|Landmanlll98ll ). A correlation is thus ex- 
pected to be present between the Call H and K index 
and the H a index, while these two indexes will in princi- 
ple be anti-correlated with the He I D3 activity indicator. 
This can be seen in Fig.[TJ though from this plot is is also 
clear that the two former lines are more sensitive to activ- 
ity variations than the latter one. This is partially caused 



by the weakness of the Hei D3 line, which it is dilluted 
by several atomic lines in the s ame region, and affected 
by telluric abso rption features (|Danks fc Lambertl 11985 ; 
ISaar et al.lll997t ). 

The wavelength region where the Can H and K lines 
are present (blue) often has a low S/N, making it difficult 
in some cases to derive a reliable activity index. Beyond 
this, and because the spectra were collected with a si- 
multaneous ThAr wavelength calibration spectrum, con- 
tamination from the nearby ThAr order may in some 
cases alter the results. To avoid this, only spectra with 
an S/N > 20 near 4000 A were used to derive the Can H 
and K "Sharps" index. 

5. Results 

5.1. Activity cycle variations 

In Figs.[5]and|3]we present time series of the three activity 
indexes for the 8 stars in our progranQ. Each point in 
the figure denotes the average of the measurements over 
one observing night (if multiple spectra existed for a given 
night). The error bars were computed using the dispersion 
of the different measurements in case several spectra per 
night were taken, or the individual error for nights when 
one single spectrum was obtained. 

A general look at the plots show that several of 
the stars present clear long-term activity cycle vari- 
ations, as expected. In general also, these variations 
are seen in all the activity indexes (see also discus- 
sions in [Livingston et al.l 120071 : ICincuneeui et al.l 12007 : 



Meunier fc Delfossell2009T . though the Call S index 



vari- 
ations seem clearer and present smaller dispersion. Our 
data suggest, however, that all these indexes can be used 
to trace long-term activity variations in solar-type stars. 

As mentioned above, the Hei index is generally anti- 
correlated with the H Q and Call H and K indexes (the 
last two are generally correlated). This is expected from 
studies of the solar active regions (|Landmanlll98ll) . The 



exception for this is the late F dwarf HD 216385. Though 
not the scope of the present paper, our results thus sug- 
gest that in F dwarfs the physics of the formation of the 
He I line is different from the one seen for early-K dwarfs. 
This star presents, however, a very stable activity level 



3 The calibration yealds Smw = 30.985 xSharps + 0.042 



4 Data available in online Tables 5, 6, and 7 
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Fig. 2. Activity indexes for HD4628, HD 16160 HD26965A, and HD 32146 as a function 
variations in the Mount Wilson H and K "S" index, middle panels the H Q index, and the 
The dashed line represents a linear fit to the data. For comparison reasons, the x-scale 
plots. 



500 1000 1500 

JD (-2453000) 

of time. Left panels denote 
right panels the He I index, 
was set constant in all the 



througout all the se ries of measurements, confirming its 
stability as shown in iBaliunas et al.l (|1995l ) - it was actu- 
ally included in our sample as a "standard" . 



the discussion of these high frequency variations to a dif- 
ferent paper. 



Some of the stars show higher frequency structure. 
This is likely due to the appearance and disappearance 
of activity phenomena such as spots and plages, and their 
modulation with stellar rotation. In the present paper we 
are mostly interested in long-term variations. We will leave 



HD219834A presents a clear spike in the H Q and He I 
indexes near JD^2453600. As shown in Sect l5.21 this star 
is a long-period binary with an eccentricity significantly 
different from zero. A fit to the radial velocities shows that 
the periastron passage occurs near the above date. We may 
thus be seeing activity induced by the fainter secondary 
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star on HD 219834A. This issue is beyond the scope of the 
present paper, but will be discussed in a separate work. 



5.2. Multiple stars 

Three of the stars in our sample are members of multiple 
stellar systems. To be able to study the influence that 
magnetic cycle variations may have on the measurement 
of RVs, we need to subtract the variations caused by the 
orbital period. In this section we describe the procedures 
used in each case. 



5.2.1. HD 16160 

HD 16160 (G1105A) is member of a known triple sys- 
tem, with G1105C orbiting G1105A with a period of 
j 60 years, orbital ecc entricity of ^0.75, and with semi- 



major axis of ~15 AU ijGolimowski et ah 2000h . Since our 



data only cover a small part of the orbital phase, we used 
the radial- velocity measurements of this star presented by 
Golimowski et al|£| , together with our own measurements, 
to fit a Keplerian function to the system and obtain a 
global orbital solution (see Tableland Fig.H]). In this pro- 



Taken directly from their postscript figure. 
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Table 4. Orbital elements of the fitted orbits for HD 16160, HD 160346, and HD 219834. See text for more details 



Parameter 


HD 16160 


HD 160346 


HD 219834 A 


Units 


P 


219001 


83.7288±0.0007 


2338±17 


[d] 


T 


2 457 280±271 


2 454 590.40±0.02 


2 453 768±19 


[d] 


a 


14.0 


0.36 


4.10 


[AU] 


e 


0.75t 


0.2048±0.0004 


0.184±0.005 




V r (HARPS) 


25.53±0.05 


21.941±0.002 


10.56±0.03 


[kms -1 ] 




144±2 


140.1±0.1 


213.9±0.7 


[degr] 


Ki 


1052±194 


5 691±3 


6017±11 


[ms- 1 ] 


o(0 - C) 


1.51 


7.11 


25.8 


[ms- 1 ] 


N 


78 


33 


25 




m2 sini 


74 


101 


448 


[M Jup ] 



f Fixed according to Golimowski ct al. (2000) 



cess we fixed the orbital period and eccentricity to the 
values mentioned above. We used the residual radial ve- 
locities of this fit for the rest of the paper. 

Leaving all the parameters free we find a slightly bet- 
ter solution with an orbital period ^11000 days, but this 
solution is not significantly better than the one adopted 



Since our data, together with the one of lGolimowski et al 



(2000), do not cover one entire orbital period, we have no 
way to constrain the orbital solution using only radial ve 
locities. 

5.2.2. HD 160346 



Halbwachs et al.l (|2003l ) used CORAVEL radial-velocity 
measurements with Hipparcos (|ESAlll997l) astrometry to 
conclude that HD 160346 (GJ 688) is a spectroscopic bi- 



nary (SB1) with a period of 83.7 days, eccentricity of 
0.20, and mass ratio M2/M1 of 0.5. A Keplerian fit to 
the HARPS radial velocities perfectly confirm this result 
(Table0] and Fig. [5]) . The residuals around this Keplerian 
fit will be used for the rest of this work. 



5.2.3. HD219834A 

HD 219834A is a known hierarchical triple system, where 
the A comp onent has been identified as a spectroscopi c 
binary SB1 (|Duauennov fc MavorlH99ll lTokovininlll997!) . 
Our HARPS data perfectly confirm (and refine) the pre- 
viously determined orbital parameters of the system. Our 
radial-velocities are best fit with a Keplerian function 
with a period of 2337 days, eccentricity 0.18, and semi- 
amplitude K^e.OlSkms" 1 (TableHand Fig.©. This sig- 
nal is compatible with a low mass (0.43 M - mini- 
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Fig. 4. Radial velocities of HD 16160 and best-fit 

Keplerian orbit. The lower panel presents the residuals Fig. 5. Phase folded radial-velocities of HD 160346 and 



of the fit. 



best-fit Keplerian orbit. 



8 



Santos et al.: Do stellar magnetic cycles influence the measurement of precise radial velocities? 




Fig. 6. Radial velocities of HD219834A and best-fit 
Keplerian orbit. The lower panel presents the residuals 
of the fit. 



mum mass) stellar companion orbiting the slightly evolved 
1.25 M star HD 219834 A. 

Unfortunately, the high residual velocities around the 
orbital fit (rms of ^25 ms -1 ), together with their struc- 
ture, strongly suggest that we are observing a blended 
spectrum and that HD 219834 A is an SB2. Alternatively, 
additional companions could be present in the system. 
Better fits are indeed obtained with 2 or 3 Keplerian func- 
tions, but because our data barely covers one orbital pe- 
riod it does not allow us to find any convincing result. 

These facts make it impossible to use the measure- 
ments of this star to study long-term and low-amplitude 
RV variations. We have thus excluded the data of 
HD 219834 A for the rest of the paper. 

5.3. HD 152391 

HD 152391 is the most active star in our sample. As ex- 
pected, a look at the obtained RVs reveal a high noise 
level close to 18ms _1 . A clear correlation between the ve- 
locities and BIS is also see n (Fig.[3 ■ As for o ther similar 
stars (such as HD 166435 - Queloz et al. 2000() . this seems 
to indicate that the observed radial- velocity variations are 
being induced by photospheric features like spots. In con- 
trast to the case of HD 166435, however, we could not find 
any clear periodicity in the data. 

We tried to correct the radial-velocities using the rela- 
tion bet ween BIS and RV, a procedure already successfully 
used bv iMelo et al.l ( 2007 1. However, the rms around the 
average velocity only decreased slightly to ~12.5ms _1 . 
This high value makes the search for any low-amplitude 



Fig. 7. Radial velocities vs. BIS for HD 152391. 

and long-term signal very difficult. We thus decided to 
keep this star out of the rest of the discussion. 

5.4. Activity, radial-velocity, and CCF parameters 

After removing HD 152391 and HD219834A from the 
list (see discussion above), we are left with 6 objects 
(HD4628, HD 16160, HD26965A, HD 32147, HD 160346, 
and HD 216385) for which we can study the influence of 
activity cycle variations on the measured radial velocities. 

We searched the literature for references to multiplicity 
among these targets. Except for HD 16160 and HD 160346 
(see discussion above), none of the remaining stars has any 
reference for a short-period (up to a few years) compan- 
ion. HD 26965A is member of a known triple system, with 
a long period of -8 000years (|Heint d [l97l . HD 216385 
is also known to be member of a common proper mo 
tion pair with the second star situated 250 arcsec away 
(jLepine fc Bongiornol 120071 ): at a distance of ^26 pc, this 
corresponds to a projected separation of ^6500 AU. For 
HD4628, HD26965A, HD 32147, and HD 216385 we thus 
have no indication that radial-velocity variations induced 
by a stellar companion can be detectable in our 5-year 
baseline. 



5.4.1. Activity and CCF parameters 

In Fig. [5] we present the yearly average time series of the 
Smw, radial- velocity, BIS, and the CCF's contrast and 
FWHIvf! In the plot, the dashed lines represent linear 
fits to the data. Each point in the plot corresponds to 
the average of the values in time bins corresponding to 



6 Data available in online Table 8 
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Fig. 8. Yearly average Smw (in each top panel), radial-velocity, BIS, and the Cross-Correlation Function's Contrast 
and FWHM function of time for HD4628, HD 16160 and HD26965A, HD 321467, HD 160346, and HDHD216385. For 
comparison reasons the x-scale was set constant in all the plots. 
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the years 2003 to 2008. The error bars are computed as 
the rms/yN, where N is the number of measurements 
in each window. The average Julian Date for each bin 
was considered for plotting reasons. The third date for 
HD 160346 (corresponding to the year 2006) has no error 
bar for FWHM and contrast, since only one measurement 
of this star was available in that year (no error bars are 
provided by the HARPS pipeline for these two parame- 
ters). For HD 32147, none of the available spectra obtained 
in 2003 was good enough for deriving a reliable "S" value, 
though radial-velocity measurements exist. 

The first remarkable result that comes out of the plots 
is that in most cases there seems to be a clear positive 
correlation between the Smw index and both the FWHM 
(the only exceptions are HD 32147 and our "standard" 
HD 216385) and the BIS of the HARPS CCF. In other 
words, the FWHM and the BIS seem to be good indica- 
tors of the stellar activity level. The opposite trend (anti- 
correlation) is seen between Smw and the CCF's contrast. 
For HD 26965 A these relations are particularly clear since 
the curvature observed in the Smw variation is also clearly 
observed in the different CCF parameters. 

The reason for these facts is not fully clear, but they 
suggest that these CCF parameters can be used as a proxy 
for the chromospheric activity level of a solar-type star. 
In particular, they can be used to follow the activity cycle 
variations of the targets in complement to measurements 
of different activity ind exes. 



Boisse et al. ( 20091 ) have shown that the parameters 



of the CCF, in particular its contrast, vary as a func- 
tion of the activity level of the star for the very active 
planet host HD 189733. The CCF appears to be shallower 
when the star is more active. These authors explain this 
variation by dark spots (typically associated with stellar 
activity) that have a spectrum that is different from the 
one emmited by the remaining stellar photosphere. These 
features will thus naturally induce changes in the overall 
spectrum that are noticeable in the shape of the CCF. We 
should note, however, that spots also influence the wings 
of spectral lines (decreasing their intensity), thus decreas- 
ing their FWHM. In this case an anti-correlation between 
FWHM and activity level could be expected. 

The resu l ts pre sented here confirm the conclusions of 
Boisse et al. ( 20091 ). Interestingly, however, our results fur- 
ther suggest that these changes are detectable even in in- 
active (solar-like) stars, such as the 6 objects discussed in 
this section. 

Variations in the CCF contrast and FWHM can also be 
induced when the cores of stronger lines start to fill as the 
star becomes more active. This will make the lines (and 
consequently the CCF) shallower, and since the FWHM is 
measured at a lower level (closer to continuum), its value 
will also increase. 

We add that a change in the shape of the CCF can 
also be induce d by variations in the convection pattern 
(|Dravinsll982t) . Observations show that solar line bisectors 
smaller lower velocity spans and convective blueshifts i n 
active regions ( Livingston 1982 ; Brandt &; Solanki 1990l ). 



This could explain the observed variations in the shape of 
the CCF as a function of the activity level of the star. 

5.4.2. Activity and radial velocity 

If on the one side activity cycle variations seem to in- 
duce clear signals on the above discussed CCF parameters 
(FWHM, Contrast, and BIS), the situation regarding the 
RV is not so clear. For HD 26965 A we find a positive cor- 
relation between the "Smw" index (that varied by —0.04) 
and the RV. Interestingly, there is a slight indication that 
a change in slope of the RV time series exists accompa- 
nying the observed variation in "S" . The same positive 
correlation is observed for HD4628, though in this case 
the RV seems to present a higher dispersion. This may 
also be caused by HD4628 presenting a lower amplitude 
variation in its activity level (A S —0.02, between the first 
and last years of measurements). 

The opposite trend, however, is found for HD 16160, 
where the "S" index increased monotonically by about 
0.054 during the period of our measurements. For this 
case we cannot fully exclude that a slightly wrong orbital 
fit could have left any trend in the residuals. 

For HD 32147 (AS -0.04) and HD 160346 
(A S —0.07), both showing clear and monotonic variation 
in the "S" index, no significant RV variation is observed 
in our data. Finally, for HD 216385, only a small "Smw" 
variation is observed (0.006). We note that this star shows 
very small long term chromospheric activity variation. It 



jmw 



3tually includec 


in our 


(Baliunas et al. 


1995). 



This is also the only F-type 
star in our sample. 

As a test, we computed the Spearman rank correla- 
tion coefficient for the "Smw" vs. RV relation. A Monte- 
Carlo test, where we generated random RV and "Smw" 
values with the same rms as the original sample shows 
that, in most cases, a higher correlation can be found in 
more than 10% of the random samplings. The only excep- 
tion is HD 26965 A, for which a false alarm probability of 
1.5% exists. The fact that we only have a few data points 
precludes any firm conclusions on this. 

In face of the present (and contradictory) results, we 
cannot conclude anything about the existence of radial- 
velocity variations induced by long-term changes in the 
chromospheric activity level in our sample of early-K 
dwarfs. More data is clearly needed, and in particular we 
need to cover the whole magnetic cycles of the target stars 
better. The study of a larger sample could also be impor- 
tant. In any case, our results suggest that long-term vari- 
ations significantly above — 1 ms _1 can be excluded in our 
targets as being caused by variations in the chromospheric 
activity level during the stellar magnetic cycle. 

We also cannot exclude the possibility that some of the 
observed RV variations are being caused by the presence of 
long-period companions, either planets or stellar in nature. 
A deep adaptive optics search for close companions to our 
targets may be crucial for clarifying this aspect. 
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6. Discussion and conclusions 

In this paper we present the results of a long-term project 
to investigate the effect that long-term stellar magnetic 
cycles may have on the measurement of precise radial ve- 
locities. For this we observed a sample of 7 late-G or early- 
K dwarfs and one late F dwarf (an activity "standard") 
for more than 5 years with the HARPS spectrograph. The 
obtained spectra allowed us to derive a precise value from 
each spectrum for the RV, a measurement of the chomo- 
spheric activity level using 3 different spectral indexes 
(Call "S", H Q , and Hei), as well as several parameters 
of the CCF (BIS, FWHM, and contrast). The results of 
this survey suggest that in our sample: 

— All the three activity indexes measured are good trac- 
ers of stellar activity level variations along the mag- 
netic cycle, though the Call "S" index presents the 
smallest dispersion. In general, the Call "S" and H Q 
indexes are correlated with each other, while the He I 
is anti-correlated with the first two. 

— The activity level variations are also clearly reflected in 
the values of BIS, FWHM, and contrast of the HARPS 
CCF. The values of BIS are particularly sensitive to 
activity level variations. This suggests that measure- 
ments of these parameters can be used to clearly diag- 
nose long-term variations in the chromospheric level of 
a star. The measurement of BIS, FWHM, and contrast 
may even likely be useful to correct the velocity mea- 
surements for the effe ct of long-term stellar activity 
( Saar fe Fischer! l2000h . in a similar way to one used 
by Melo et al.1 ( 20071 ) to correct for higher frequency 
noise. These results may be of extreme importance for 
present and future high-precision RV planet searches. 

— Although some of our targets show hints of low- 
amplitude (at the ~lms _1 level) and long-term RV 
variations that could be caused by variations in the 
activity level of the star, our results are not conclu- 
sive about the nature and amplitude of this effect. Our 
data suggests that for early-K dwarfs the variation in 
the stellar activity level along the magnetic cycles does 
not strongly induce variations in the measured RVs. 
Early-K dwarfs are thus bona fide targets to search for 
very low-mass planets using precise radial velocity in- 
struments. We note that these targets already have the 
lowest granulation and oscillation "noise" level among 
solar type stars (Dumusque et al., in preparation). 

A better understanding of the influence of long-term 
stellar magnetic cycles on the radial- velocities is certainly 
crucial to plan future radial- velocity planet searches with 
a new generation of instruments like ESPRESSO@VLT0 
or CODEX@E-ELT (jPasauini et al.ll2008l) . These instru- 
ments will be able to obtain RVs with a precision of better 
than 10 cms -1 , allowing the detection of Earth-like plan- 
ets in the habitable zone of nearly solar-type stars. 

To further investigate the above discussed issues, 
a continuation of the present program is needed. 

7 http://espresso.astro.up.pt 



Furthermore, it would be very important to extend our 
sample to earlier type G-dwarfs (more similar to our Sun) . 
Convectiv e blueshifts are lower in K-dwarfs than in F- and 
G-dwarfs (|Gravlll992f) . In active regions the convective ve- 
locities are lower , also implying smaller li ne asymmetries 
(|Livingstonl 119821 iBrandt fc Solankilll990t) . We may thus 
expect a stronger eff ect of chrom ospheric activity for ear- 
lier type dwarfs (e.g. Saar 20091 ). 

As a complement, it would also be important to inves- 
tigate the existence of long period stellar companions in 
detail, which could be able to induce long-term RV varia- 
tions. An adaptive optics survey should thus be conducted 
to investigate this possibility. 
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